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Abstract 

The exclusive rare radiative -B-decays are studied in the framework of the relativistic quark model based 
on the quasipotential approach in quantum field theory. The large recoil momentum of the final vector 
meson allows for the expansion of the decay form factor in inverse powers of the b-quark mass. This 
considerably simplifies the analysis of these decays. The l/rrifc expansion is carried out up to the second 
order. The form factor of B — > K*j decay is found to be F-y ~* 7 (0) = 0.32 ± 0.03 and it leads to 
BR(B — > K*j) — (4.5 ± 1.5) x 10 -5 , which is in agreement with the recent CLEO data. The form factors 
of the decays B — > ,07 , B s — ► cf>j and B s — > K*^ are also considered. The relation between rare radiative 
and semileptonic i?-decays into light vector meson is discussed. 
PACS number(s): 12.39Ki, 13.20He, 13.40Hq 



1 Introduction 

Rare radiative decays of B-mesons represent an important test of the standard model of electroweak 
interactions. These transitions are induced by flavour changing neutral currents and thus they are sen- 
sitive probes of new physics (see e. g. [1]). Such decays are governed by one-loop (penguin) diagrams 
with the main contribution from virtual top quark and W boson. Therefore, they provide valuable in- 
formation about the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements Vtd, Vts and Vtb, and about 
the top quark mass. The first exclusive decay B — * K*j has been observed by CLEO [2]. The mea- 
sured branching ratio is BR(B — > K*j) = (4.5 ± 1.5 ± 0.9) x 10~ 5 . Also the inclusive decay rate 
B — ► A s 7, obtained from the measurement of the photon energy spectrum, has been reported by CLEO: 
BR(B -> X sl ) = (2.32 ± 0.51 ± 0.32 ± 0.20) x 10" 4 [3]. 

Theoretical analysis of rare radiative decays is based on the effective Hamiltonian, which is obtained 
by integrating out the heavy particles. The renormalization of the effective Hamiltonian coefficients has 
been calculated to leading order [4,5]. It turns out that renormalization effects play an important role. 
They increase the decay amplitude approximately by a factor of two [4]. Some of the next-to-leading 
order corrections are also known [6,7]. The hadronic matrix elements of the effective Hamiltonian for 
inclusive rare radiative decays have been calculated in the framework of the heavy quark effective theory 
(HQET) [8] . It has been shown that the leading term in 1/ mj expansion corresponds to the parton model 
prediction and the nonperturbative corrections are suppressed by 1/m^ factor [9,10]. For the calculation 
of the exclusive decay rates it is necessary to know the relevant hadronic form factors. In the case of 
B — > AT*7 decay only one hadronic form factor -Fi(O) contributes. It has been calculated in the framework 
of nonrelativistic quark model [11] and QCD sum rules [12-14]. In [15] the heavy quark mass limit has 
been applied to exclusive B — > K*j decay. However, the s-quark in the final K* meson is not heavy. Its 
mass is of order of A parameter, which determines the scale of 1/mq corrections in HQET [8]. Thus very 
substantial corrections to this limit come from the whole series in l/m s . Nevertheless, the ideas of heavy 
quark expansion can be applied to the exclusive B — > K*^ decay. From kinematical analysis it follows 
that the final K* meson bears large relativistic recoil momentum |A| of order of rrib/2 and the energy 
of the same order. So it is possible to expand the matrix element of the effective Hamiltonian both in 
inverse powers of b-quark mass from the initial B meson and in inverse powers of the recoil momentum 
I A I of the final K* meson. As a result the expansion in powers of 1/nib arises. The aim of this paper is 



to realize such expansion in the framework of the relativistic quark model. We show that this expansion 
considerably simplifies the analysis of exclusive rare radiative decays of B mesons and reduce the model 
dependence of the result. The leading, first and second order terms of the l/mj expansion are calculated. 
It is necessary to note that rare radiative decays of B mesons require the completely relativistic treatment, 
because the recoil momentum of final vector meson is highly relativistic. 

Our relativistic quark model is based on the quasipotential approach in quantum field theory with 
the specific choice of the qq potential. It provides a consistent scheme for calculation of all relativistic 
corrections at a given order of v 2 /c 2 and allows for the heavy quark I/tuq expansion. This model has been 
applied for the calculations of meson mass spectra [16], radiative decay widths [17], pseudoscalar decay 
constants [18], scmileptonic [19] and nonleptonic [20] decay rates. The heavy quark I/tuq expansion in 
our model for the heavy-to-heavy semileptonic transitions has been developed in [21] up to l/rrig order. 
The results are in agreement with the model independent predictions of HQET [8]. The rare decay 
B — > K*j has been considered in our model in [22]. Here we refine our previous analysis with more 
complete account of relativistic effects and using l/m& expansion. We also consider some other exclusive 
radiative decays, including the decay B s — > (j)-f and the CKM-suppressed decays B pq and B s — > K*j. 

The paper is organized as follows. In Sect. 2 we briefly review the theoretical predictions for the 
inclusive rare radiative B-decays and define the form factor which governs the exclusive decays. The 
relativistic quark model is described in Sect. 3, and in Sect. 4 it is applied for the calculation of the rare 
radiative decay form factor F±. The 1/rrib expansion for this form factor is carried out in Sect. 5. Our 
numerical results for the form factors of the decays B — * K*j , B s — > <fry , B — > p'y and B s — > K*-f are 
presented in Sect. 6. We also discuss the relations between rare radiative and semileptonic _B-decays into 
light vector mesons. Sect. 7 contains our conclusions. 



2 Inclusive and exclusive rare radiative 5-decays 

In the standard model, B-decays are described by the effective Hamiltonian, obtained by integrating out 
the top quark and W boson and using the Wilson expansion [4] . For the case of b — > s transition: 

AQ 8 

H ef f(b ^a) = -^-V t * s V tb J2CMOM, (1) 

where Vij are the corresponding CKM matrix elements, {Oj} are a complete set of renormalized dimension 
six operators involving light fields, which govern b — > s transitions. They consist of six four-quark 
operators Oj (j = 1, . . . , 6), which determine the non-leptonic -B-decay rates, the electromagnetic dipolc 
operator O7 

7 = j^sa^(m b P R + m s P L )bF l _ lv , P R , L = (1 ± 75 )/2, (2) 

and the chromomagnetic dipole operator 0$, which are responsible for the rare -B-decays B — > X s j and 
b — > s + g, respectively [4]. The Wilson coefficients Cj(/i) are evaluated perturbatively at the W scale and 
then they are evolved down to the renormalization scale fi ~ m& by the renormalization group equations. 
The coefficient of the magnetic operator, C7, has been calculated to leading logarithmic order [4,5]. The 
next-to-leading order corrections to the anomalous dimension matrix are also partially known [6]. 

The dominant contribution to the inclusive decay width T(B — > X s -f) comes from the magnetic 
moment term CV(/z)07(/x). Thus it is necessary to calculate the hadronic matrix element of this operator. 
Recently it has been observed that the matrix elements of inclusive decays of hadrons containing a heavy 
quark Q allow for a systematic expansion in powers of 1 /mg [9] . The leading order term of this expansion 
reproduces the parton model rate and the nonperturbative corrections appear only at the second order 
of 1/rriQ expansion [9]. Therefore, the decay rate for B — > X s -f is [10] 
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with p — (3mb + 5m s ) / (3rrib — 3m s ). The parameter Ai is related to the kinetic energy of the 6-quark 
inside B meson and A 2 is related to the B — B* mass splitting [8]. This rate is strongly dependent on 
the value of 6-quark mass mb, which depends on the definition. So it is more convinient to connect 
T(B — > X s -f) with the experimentally observed semilcptonic decay rate. 

The branching ratio BR(B — > X s j) can be expressed in terms of the inclusive semileptonic branching 
ratio BR(B -» Xlv t ) as [7] 

BR(B - X sl ) = 6 2 'Y'f - |gr(r» t )Wr» t ) ^ 

where g(r) is the phase-space factor for Tip — > c£^): <?(r) = 1 — 8r 2 +8r 6 — r 8 — 24r 4 ln(r), and experimental 
value BRiB — » Xtvt) k, 10.5% is used for the semileptonic branching fraction. The function /(r) accounts 
for QCD corrections to the semileptonic decay [23] and, for a typical quark mass ratio of r = 0.35 ± 0.05, 
it has the value f(r) — 2.37 =F 0.13. The factor Kim) contains 0(a 3 ) corrections to the B — ► X s j rate, 
due to the QCD bremsstrahlung processes [7]. The resulting branching ratio is [7,24] 

BRiB — > X s j) = (3.0 ± 1.2) x 10~ 4 , (5) 

for m t in the range 100 < m t < 200 GeV. 

This result agrees with a recent data on the inclusive decay rate by CLEO collaboration [3]: 

BRiB -> X s -f) = (2.32 ± 0.51 ± 0.32 ± 0.20) x 10~ 4 . (6) 

In what follows we shall discuss the exclusive rare radiative decays of B mesons, like B u ^ — > K*j, 
B s — » 07 and the CKM-suppressed chanels -B Ul d — > P7, £> s — > if* 7 . We denote all these modes by 
B — ► V7 , where i? is or _B S meson and V is if*, cj> or p. The main contribution to the exclusive 
B — > K7 decay amplitude comes from the magnetic moment operator O7. The relevant matrix element 
has the covariant decomposition [11]: 

(Vi Pv ,e)\fia^P R b\Bi PB )} = ie^e^p^F^q 2 ) 

+ [e;(Mi - My) - (e* • q)i PB +Pv),]G 2 iq 2 ), (7) 

where G2(0) — Fi(0)/2, e is a polarization vector of final vector meson, q = pb — Pv denotes the 
four-momenta of the emitted photon. The exclusive decay rate is given by 

TiB -» if* 7 ) = ^\ V ; s V tb \ 2 \F im 2 \C 7 irn b )\ 2 i m t + m 2 ) (M| M f K ' ] \ (8) 

and the analogous expressions can be written for T(i? s — > ^7), T(_B — > ^7) and T(i3 s — > K*-f). Then the 
exclusive branching ratio BRiB — ► if *7) is 

BRiB ^ K*~f) = R- BR(B ^ X s j), (9) 

where the ratio of the exclusive to inclusive radiative decay rates is determined by (3), (8) to be: 

R = E(g ^ jgy) _ (1 - M|,/M 2 ) 3 M| |fi(0)| 2 f , 

r(B^X s7 ) (l- m 2/ m 2) 3 if(m,)(l + (A 1 -9pA 2 )/(2m 2 ))- 1 ! 

The factor if(m 6 ) w 0.83 [7], the parameters Ai = -0.30 ± 0.30 GeV 2 [8,25] and A 2 = 0.12 ± 0.01 GeV 2 
[8]. 

We shall calculate the form factor F\ (0) with the account of l/m& corrections up to the second order 
in the framework of the relativistic quark model based on the quasipotential method. 
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3 Relativist ic quark model 



In the quasipotcntial approach meson is described by the wave function of the bound quark-antiquark 
state, which satisfies the quasipotential equation [26] of the Schrodinger type [27]: 

TS? - Q * M(P) " / HrnP,q;M)*„( q ), (ii) 

where the relativistic reduced mass is 

_ M*-{ml-ml)\ 
m ~ AM* ' ^ 

6 2 (M) = [M 2 - {ma + m b ) 2 ][M 2 - (m a - m b f] ^ ^ 

m a ,b are the quark masses; M is the meson mass; p is the relative momentum of quarks. While construct- 
ing the kernel of this equation V(p, q; M) — the quasipotential of quark-antiquark interaction — we have 
assumed that effective interaction is the sum of the one-gluon exchange term with the mixture of long- 
range vector and scalar linear confining potentials. We have also assumed that at large distances quarks 
acquire universal nonperturbative anomalous chromomagnetic moments and thus the vector long-range 
potential contains the Pauli interaction. The quasipotential is defined by [16]: 

V{p,q,M) = u a {p)u b {-P){^ S D^{k)^ + V c v onf {k)T^ + V c s onf {k)}u a {q)u b {-q), (14) 

where as is the QCD coupling constant, is the gluon propagator; 7^ and u{p) are the Dirac matrices 
and spinors; k = p q; the effective long-range vector vertex is 

r„(k) = 7m + ^/*^, (15) 

k is anomalous chromomagnetic quark moment. Vector and scalar confining potentials in the nonrela- 
tivistic limit reduce to 

V c v onf {r) = {l-e){Ar + B), Vg nf {r) = e{Ar + B), (16) 

reproducing V^°"/ eZ (r) = V^ on j + V^ on j = Ar+B, where e is the mixing coefficient. The explicit expression 
for the quasipotential with the account of the relativistic corrections of order v 2 /c 2 can be found in ref. 
[16]. All the parameters of our model: quark masses, parameters of linear confining potential A and 
B, mixing coefficient e and anomalous chromomagnetic quark moment n were fixed from the analysis 
of meson masses [16] and radiative decays [17]. Quark masses: m b — 4.88 GeV; m c — 1.55 GeV; 
m s = 0.50 GeV; m u4 = 0.33 GeV and parameters of linear potential: A = 0.18 GeV 2 ; B = -0.30 GeV 
have standard values for quark models. The value of mixing coefficient of vector and scalar confining 
potentials e = —0.9 has been primarily chosen from the consideration of meson radiative decays, which 
are very sensitive to the Lorentz-structure of the confining potential: the resulting leading relativistic 
corrections coming from vector and scalar potentials have opposite signs for the radiative Ml - decays 
[17]. Universal anomalous chromomagnetic moment of quark n = — 1 has been fixed from the analysis of 
the fine splitting of heavy quarkonia 3 Pj- states [16]. 

Recently we have considered the expansion of the matrix elements of weak heavy quark currents 
between pseudoscalar and vector meson states up to the second order in inverse powers of the heavy 
quark masses [21]. It has been found that the general structure of leading, subleading and second order 
1/toq corrections in our relativistic model is in accord with the predictions of HQET. The heavy quark 
symmetry and QCD impose rigid constraints on the parameters of the long-range potential of our model. 
The analysis of the first order corrections [21] allowed to fix the value of effective long-range anomalous 
chromomagnetic moment of quarks k = — 1, which coincides with the result, obtained from the mass 
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spectra [16]. The mixing parameter of vector and scalar confining potentials has been found from the 
comparison of the second order corrections to be e = — 1. This value is very close to the previous one 
e = —0.9 determined from radiative decays of mesons [17]. Therefore, we have got QCD and heavy quark 
symmetry motivation for the choice of the main parameters of our model. The found values of e and 
k imply that confining quark-antiquark potential has predominantly Lorentz-vector structure, while the 
scalar potential is anticonfining and helps to reproduce the initial nonrelativistic potential. 



4 Calculation of the rare radiative decay form factor 



The amplitude of the exclusive B — > Vj decay is proportional to the hadronic matrix element (7) of the 
magnetic moment operator O7. Thus it is necessary to calculate the transition form factor -Fi(O). 

The matrix element of the local current J between bound states in the quasipotential method has the 
form [28]: 

d 3 p d 3 q 

W) 

where T M (p, q) is the two-particle vertex function and ^v,b are the meson wave functions projected onto 
the positive energy states of quarks. 

In the case of rare radiative decays J M = j 'io '^ v k v Pub and in order to calculate its matrix element 
between meson states it is necessary to consider the contributions to T from Figs. 1 and 2. Thus the 
vertex functions look like 



(V\J,(0)\B) = j ^E^* v (p)r M (p,q)* B (q), (17) 



and 



rW(p,q)=S / (p 1 )^^^(l+ 7 5 K(g 1 )(27r)^(p 2 -q 2 ), (18) 



rf(p,q) - u s { Pl )u q { P2 ) l -Ua 1 , v (l+ll) p (tVjl 1 ] T 7iV(P2~q 2 ) 

2 L e b [ki) + e b {pi) 

+y(p 2 -q 2 ) ^{jl 1 ] r 7iW(l +7i 5 )Wgi)V<Z2), (19) 
Sf ( fc l )+£/(<?! ) J 



where ki = pi - A; k^ = qi + A; A = p B -p y ; e{p) = (m 2 + p 2 ) 1/2 ; 

e(p) - (to7° +7°(7P)) 



A<->(p) = 



2e(p) 



Note that the contribution r^ 2 ^ is the consequence of the projection onto the positive-energy states. The 
form of the relativistic corrections resulting from the vertex function T^ 2 ' is explicitly dependent on the 
Lorentz-structure of gg-interaction. 

It is convenient to consider the decay B — > V-y in the B meson rest frame. Then the recoil momentum 
A of vector meson V is highly relativistic (| A| « M B /2). The wave function of the meson moving with 
the momentum A is connected with the wave function at rest by the transformation [28] 

*va(p) = <«a)^ /2 (0%o( P ), (20) 

where D 1 / 2 {R) is the well-known rotation matrix and R w is the Wigner rotation. 

The meson wave functions in the rest frame have been calculated by numerical solution of the quasipo- 
tential equation (11) [29]. However, it is more convenient to use analytical expressions for meson wave 
functions. The examination of numerical results for the ground state wave functions of mesons containing 
at least one light quark has shown that they can be well approximated by the Gaussian functions 

* m (p) = *mo(p) = (^-) 3/4 ^(-^r), (21) 
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with the deviation less than 5%. 
The parameters are 



Pb = 0.41 GcV; fi K * = 0.33 GeV; (i p = 0.31 GeV; 
(3 Bs = 0.46 GeV; /3 = 0.36 GeV. 



(22) 



Substituting the vertex functions (18), (19), with the account of wave function transformation of 
vector meson (20), in the matrix element (17) we get for the form factor F\ the following expression: 



Fi(0) 
ff>(0) 



i^f (2) (0) 



*f (2) (o) 



F^>(0) + eif (2) (0) + (1 - e)Ff w (0), 

2e a 



(23) 



E v f d 3 p 



M B J (2tt) 3 



*v P 



E v + My ' 



lsf(p + A)+mf e b (p)+m b 



x< 1 + 



Vl + (PA) 



(e/(p + A) + m f )(e b (p) + m b ) 



V 2 £/ (p + A) 
(M B - £7 V ) 



2e 6 (p) 



£/(p + A) +m/ 



+ 



(PA) 



+ 



£/(p + A)+m/ £ 6 (p)+m 6 



pI+pI 



1 



1 



+ 



2(£y + My) \e q (p) + m q \e b (p) + m h e f (p + A) + m f 
4 



(e/(p) +m f )(e b (p) +m b ) 



>*b(p), 



(24) 



2£„ 



^l^v ^P^ 

M B 7 (2tt) 3 y V £y + My 



, £/(p + A)+m / 



f £/(A) 



I 2£ / (A)( £/ (A)+™ / ) 
2£ C 



1 + 



/ 2 £/ (p + A) 

Mb — Ey ( P A ) S 
2m h A s 



x My - £/ p + 



A - e q (p 



2e n 



Ey + My 

M B + My - e b (p) - e q (p) -e f (p + 



E v + M v 



(Mb - E v ) 



(PA) 
A 2 



2e„ 



A )-^(p + ^k7 A ) 



+ 



l 2m b (e b (A) + m b ) 

M B - My - e b (p) - e q {p) + e f (p + e7Tm7 A ) + £ i (p 



2e„ 



E v +M v 



A 



) 



Ey f d 3 p 



*V P + 



2e/(A)(£/(A)+m/ 
2£„ 



>*b( P )(25) 



My J (2tt) 3 

[ £/(A)-m / 
\2£/(A)(£ / (A)+to / ) 

2e,_ 



Ey + My 
1 + 



Sf(p + A) + nif 
^ f (p + A) 



X My - £/ [p + 



Mb — Ey (pA) \ 
2m 6 A 2 / 



2£„ 



£y + My 



£y + My 



A + (Mb - Ey) 



(PA) 



'My - e, (p + 2e q E v + My A) -e q (p+ y^-A 



A 2 

M B - e b ( P ) - e q ( P ) 



1 + K 



2m b (e b (A) + m b ) 
1 1 



2(e q (p)+m q ) \e b (A)+m b e f (A) 



2 £/ (A)(£ / (A)+m / ) 
Mb - My - e 6 (p) 
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(26) 

where the superscripts "(1)" and "(2)" correspond to Figs. 1 and 2, S and V — to the scalar and vector 
potentials of gg-interaction; 

M% — My M% + M v , 

|A| = 2M B ' V = 2M B ' ( } 

and z-axis is chosen in the direction of A. The contributions to the form factor coming from Fig. 2 
are proportional to the binding energy of B or V mesons. Taking into account that in our model 
m b + m d ~ M b ; m b + m s ~ M Bs ; m s + m d ~ M K * \ 2m s ~ M^; 2m d ~ M p and |A| 2 ~ m 2 /4 » (p 2 ) we 
have neglected the terms proportional to the product of binding energies and ratios p 2 /e^ (A), p 2 /e^(A). 
The omitted terms are of order l/m|. 

The expressions (23)-(26) for the form factor i*i(0) differ from our previous results [22] in the argument 
of the wave function of final vector meson V. In [22] we have used for simplicity the nonrelativistic limit 
p + -j^- A for this argument, taking into account the approximation of the wave functions by Gaussians. 
However, our recent analysis of the 1/mg corrections to the semileptonic decays [21] has shown, that such 
approximation is not completely adequate. It is necessary to take relativistic expression p + E ^ 9 Mv A 
in the argument of the wave function, what is done in eqs.(23)-(26). 



5 I/rut expansion for the form factor Fi(0) 

At present considerable theoretical interest is attracted to HQET. This theory is based on QCD and 
heavy quark expansion. Additional heavy quark symmetries (for a review see e. g. [8]) arise in the limit 
of infinitely heavy quark mass. HQET provides a systematic expansion in inverse powers of the heavy 
quark mass of hadronic matrix elements between mesons with one heavy and one light quarks. In the 
case of exclusive heavy-to-heavy transitions heavy quark symmetries considerably reduce the number of 
independent form factors in each order of 1/toq expansion. This simplifies the theoretical analysis of such 
decays. In exclusive heavy-to-light meson transitions (such as the B — > V7 decays) the predictive power 
of HQET is essentialy less. Only the relations between various form factors of the semileptonic and rare 
-B-decays are imposed in the limit of infinitely heavy 6-quark mass [30]. However, as it has been already 
mentioned in the introduction, the large value of the recoil momentum of final vector meson | A| ~ m b /2 
allows the expansion in powers of 1/mf, for the radiative rare decay form factor i*i(0). This expansion 
leads to considerable simplification of the formulae for i*i(0). It has been already partly used by us in 
derivation of (23)-(26). The large value of recoil momentum |A| allowed us to neglect p 2 in comparison 

with A 2 in the quark energy £f(p + A) in final meson in the expressions for f^ 2 ' and F^ 2 \ Thus 
we were able to perform one of the integrations in the current matrix element (17) using quasipotcntial 
equation as in the case of the heavy final meson [17,19,28]. As a result, we get more compact formulae. 
In this section for the sake of consistency we carry out the complete expansion of (23)-(26) in inverse 
powers of 6-quark mass. 

In HQET the mass of B meson has the following expansion in l/m b [8] 

where parameter A is the difference between the meson and quark masses in the limit of infinitely heavy 
quark mass. In our model A is equal to the mean value of light quark energy inside the heavy meson 
A = {s q )B ~ 0.54 GeV [21]. Am B arises from the first-order power corrections to the HQET Lagrangian 
and has the form [8]: 

Am| = -Ai-3A 2 . (29) 
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The parameter Ai results from the mass shift due to the kinetic operator, while A2 parameterizes the 
chromomagnctic interaction [8]. The value of spin-symmetry breaking parameter A2 is related to the 
vector-pseudoscalar mass splitting 

A 2 w J(M|, - M%) = 0.12 ± 0.01 GeV 2 . 

The parameter Ai is not directly connected with observable quantities. Theoretical predictions for it vary 
in a wide range: Ai = -0.30 ± 0.30 GeV 2 [8,25]. 

In the limit wlq — > 00, meson wave functions become independent of the flavour of heavy quark. Thus 
the Gaussian parameter (3b in (21) should have the following expansion [21] 

Pb = P - — + O (^-] , « 0.42 GeV, (30) 

where the second term breaks the flavour symmetry and in our model is equal to A/3 2 w 0.045 GeV 2 [21]. 

Substituting (28) in (27) we get the l/m b expansion of the recoil momentum and the energy of final 
vector meson: 

^-M*]] + o(' 



|A| 


m b 
2 


(- 


m b 


1 

m 2 


Ey 


m b 
2 


(- 


m b 


1 

to 2 



Now we use the Gaussian approximation for the wave functions (21). Then shifting the integration 
variable p in (23)-(26) by — g + M A, we can factor out the A dependence of the meson wave function 
overlap in form factor F± . The result can be written in the form 

F 1 (0)= ^i(A 2 )cxp(-CA 2 ), (32) 

where | A| is given by (27) and 

2 

[Pl+Pl) \M B +M V/ 
here A is equal to the mean value of light quark energy between heavy and light meson states: 

A = <£,). (34) 

Expanding (33) in powers of l/m b we get 



ca2 _ 2A 2 A 2 (/? 2 + [3 v )(E v + M V ) 2 2A 2 ^ M B - M v ^ ^ 



9 A 2 / My\ „ ( 1 



CA 2 = — r, 1 - 4^ + O — , (35) 

p 2 V m b / \ m b J 

2d 2 

where 77 = ^ ^ ^ . We see that the first term in this expansion is large. Really, even in the case of the 

lightest final meson p: 4M p /m b w 0.63. The value of this correction to the form factor Fi(0) is also 
increased by the exponentiating in (32). Therefore, we conclude that the first order correction in l/m b 
expansion of i*i(0), arising from the meson wave function overlap, is large. Thus, in the following, we 
use uncxpanded expression (33) in the exponential of the form factor fi(0) in (32). 

In contrast to the meson wave function overlap the factor ^(A 2 ) in (32) has a well defined l/m b 
expansion. First and second order corrections are small. Substituting the Gaussian wave functions (21) 
in the expressions for the form factor (23)-(26), with the value of anomalous chromomagnetic quark 
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moment k = — 1, and using (32) and the expansions (28)-(31), we get up to the second order in l/m& 
expansion: 

^(A 2 ) = J' 1 (1) (A 2 )+e^f (2) (A 2 ) + (l-e)^ (2) (A 2 ); (36) 



2wfc V \ \ 2 + m q 3 £/ + to/ 
+i( M v "•»? - f<P 2 > - I A2 " 2 + A "< 4 "'/ + 2M >'> 

4 (iTrW)* 2 ^ -*»'»)) : (37) 



-5(2)/ - ..1 / 1 1 t 1 r 



^^(A 2 ) = TV x — 1- — (3m / + A+ -At,)) (M V -(£,) -(£/» 

i 



f-— t=(A + M V - (e f ) - 2<e,» + 2(A -M V + (e f )) 
3 /p 2 \ 1/5 1 \ /p 2 



3 V2 2 + VEJ \e/, 
^ (2) (A 2 ) = tf±(-L (l--L(3m, + A + lA,)) (M v -(e- g ) - <e,» 



.h(. 



i 



= (M v -<e / )-<e,))+2(A -<e,)) 



m 2 \2 + V5 

If 2 \ / p 2 \ / 1 1 \ / p 2 



6 \ 2 + s/Z)\e } / \2 3(2 + V5) 



(38) 



(39) 



where N = { ^ B ^ ^ = (f^ 7 ?) i s due to the normalization of Gaussian wave functions in (21); 



£ i = yp 2 +m 2 + A 2 i] 2 (i = f,q), i. e. the energies of light quarks in final vector meson acquire 
additional contribution from the recoil momentum. The averaging is taken over the Gaussian wave 
functions of B and V mesons, so it can be carried out analytically. For example, 

&) = ■±= 7 ?$ 7 =e'K 1 (z), (40) 

where m 2 = mf + k 2 rf and K\(z) is the modified Bessel function; z = fh 2 /2rjf3y. Analogous expressions 
can be obtained for the other matrix elements in (37)-(39). 

Our final result for the rare radiative decay B — ► form factor -Fi(O) is given up to the second order 
of l/m b expansion by eqs. (32), (33) and (36)-(39). 



6 Results and discussion 

Using the parameters (22) of the wave functions (21) in the expressions for the form factor i*i(0) (32), (33) 
and (36)-(39) we get (for the values of the anomalous chromomagnetic quark moment k = — 1 and the 
mixing parameter of vector and scalar confining potentials s = — 1 [21]): 

F^ K ' J (0) = 0.32 ±0.03 F^ PJ (0) = 0.26 ± 0.03 

Fj Ba_> * 7 (0) = 0.27 ±0.03 F 1 Bs_> * r * 7 (0) = 0.23 ± 0.02. (41) 
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The theoretical uncertainty in (41) results mostly from the approximation of the wave functions by 
Gaussians (21) and does not exceed 10% of form factor values. The contributions of higher order terms 
in I /nib expansion in (36)-(39) are negligibly small. Thus the unexpanded (23)-(26) and expanded (41) 
values of form factors differ unessentially. This conclusion is confirmed by numerical analysis. 

Our results (41) for the form factor values are in good agreement with recent calculations within the 
light-cone QCD sum rule [13] and hybrid sum rule [14] approaches. The comparison of predictions is 
given in Table 1. All results agree within errors. 

The calculated value of form factor F^ K 7 in (41) yields for the ratio (10) of exclusive to inclusive 
radiative decay widths: 

T(B — > K*~/) 

R(B K*<y) = = (15 ± 3)%. (42) 

Combining this result with the QCD-improved inclusive radiative branching ratio BR(B — > X s "f) given 
by (5), we get 

BR th (B -> K*i) = (4.5 ± 1.5) x 10~ 5 . (43) 
This branching ratio agrees well with experimental measurement by CLEO [2] 

BR exp (B -> K*j) = (4.5 ± 1.5 ± 0.9) x 10~ 5 . 

Our numerical analysis has shown that 1/rrib corrections in (37)-(39) give rather small contributions 
to the decay form factor. So in the first crude approximation we can neglect them. This corresponds to 
the limit nib — ► 00 f° r -^(A 2 ). Then the dependence on the Lorentz structure of confining potential is 
lost and we get a simple formula 

with w — ttS^jt^- We have introduced the function 

which in the limit of infinitely heavy final quark in V meson, coinsides with the Isgur-Wise function of 
our model [21]: 

„ , ( 2 \ 1/2 ( \ 2 w-\\ , M 2 +M 2 -q 2 

£H = 7 ex P ^777 j w = vv'= z „ - Y r — -• (46) 

sv ; \w + lj \ /3 2 w + lJ 2M B M V 

Thus if we take the formal limit of infinitely heavy s-quark, as it is done in [15], we get for B — ► K**y form 
factor 

„ , n . M B +M K , 

^ ' = 2v ^bM k . for to 6 ^ oo, m s ^ oo, (47) 

with w = 2M B M K t ■ ^ ne re l a tion (47) is, certainly, not adequate, because s-quark is not heavy. It 
overestimates the form factor approximately by a factor of 1.5. 

The approximate formula (44) gives the values of form factors, which are slightly higher than (41). 
However, the difference does not exceed 10%. Thus we can conclude that the (^-dependence of form 
factor F\ near q 2 = is determined by the function (45). 

We can use our results for F\ to test the HQET relations [30] between the form factors of rare radiative 
and semileptonic decays of B mesons. Isgur and Wise [30] have shown that in the limit of infinitely heavy 
6-quark mass an exact relation connects the form factor F\ with the semileptonic decay form factors 
defined by: 

(V( Pv ,e)\f^(l- l5 )b\B( PB )) = -(M B + M v )A 1 (q 2 )e; 



10 



where the ellipses denote terms proportional to (ps + Pv)n or Qn- This relation is valid for q 2 values 
sufficiently close to q^ax = {Mb — My) 2 and reads: 

_ q 2 + M 2 -M 2 V(q 2 ) M B + M V , 

Fliq } " 2M B M B + M V + 2Mb Mq Y (49) 

It has been argued in [31,32,14], that in these processes the soft contributions dominate over the hard 
perturbative ones, and thus the Isgur-Wise relation (49) could be extended to the whole range of q 2 . 

We can calculate the form factors of semileptonic B-decay in vector light meson at the q 2 = point 
using the method developed here for rare radiative decay form factor Fi(0). It is easy to see that the 1/rrib 
corrections coming from the overlap of meson wave functions are also large. The similar exponential factor 
as in eq. (32) arises, and we shall not expand it in the I/771&. The l/m& corrections to the preexponcntial 
factor are again rather small. So here we limit our analysis only to the leading order of 1 jrnb expansion 
for this factor. In this approximation the semileptonic decay form factors are independent of the Lorentz- 
properties of the confining gg-potential and are given [33] by a simple formulae 

where £{w) is defined by (45). In the limit of infinitely heavy b- and /-quarks the form factors (50) and 
(51) satisfy all the relations imposed by HQET [8]. The ^-dependence of the form factors near q 2 = 
is determined for V\ by the function £{w) and for A\ by the product |(1 + w)£{w). 

It is now easy to check that form factors (50), (51) and (32) satisfy Isgur-Wise relation (49) for the 
q 2 values near zero. The complete analysis of 1/toj corrections to the semileptonic 5-decays into light 
mesons up to the second order terms will be given elsewhere [33] . 



7 Conclusions 

We have investigated the rare radiative decays of B mesons. The large value (~ mb/2) of the recoil 
momentum of the final vector meson requires the completely rclativistic treatment of these decays. On 
the other hand, the presence of large recoil momentum in the energy of the final meson allows for the 
1/mfe expansion of weak decay matrix element. The contributions to this expansion come both from 
heavy &-quark mass and large recoil momentum of light vector meson. 

Using the quasipotential approach in quantum field theory and the relativistic quark model, we 
have performed the 1/mj expansion of the rare radiative B-decay form factor i*i(0) up to the second 
order. We have found that 1/rrib corrections coming from the overlap of meson wave functions arc 
large. This agrees with the results of the hybrid sum rules [14], where considerable l/m& corrections 
for the form factor -Fi(O) have been also obtained. Thus we treat the A 2 -dependence of F±, arising 
from the overlap of meson wave functions, without expansion. All other l/m& corrections turn out to 
be small. They reduce the value of the form factor i'i(O) approximately by 10%. We have calculated 
the form factors for the decays B — ► K*"f , B — > pj and B s — > (jyy , B s — * K*^ . The found values 
(41) are in good agreement with the predictions of QCD sum rules [13,14]. Combining our result for 
Ff^ K 7 (0) = 0.32 ± 0.03 with the QCD-improved theoretical prediction for inclusive branching ratio 
(5), we obtain BR(B — ► K*~f) = (4.5 ± 1.5) x 10 -5 , which is in accord with the experimental branching 
ratio BR(B -> K*-y) = (4.5 ± 1.5 ± 0.9) x 10~ 5 [2]. 

We have also evaluated [33] the form factors of B meson semileptonic decays into light vector mesons 
at the point of maximal recoil. It has been found that the relations between rare and semileptonic decay 
form factors [30], obtained in the limit of infinitely heavy 6-quark, are satisfied in our model. 
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Table caption 

Table 1 Comparison of our predictions for the rare radiative decay form factor i*i(0) with light-cone 
QCD sum rule [13] and hybrid sum rule [14] results. 

Figure captions 

Fig. 1 Lowest order vertex function 

Fig. 2 Vertex function with the account of the quark interaction. Dashed line corresponds to the effective 
potential (14). Bold line denotes the negative-energy part of quark propagator. 



Table 1 



Decay 


our results 


[13] 


[14] 


B -> K*-y 


0.32 ±0.03 


0.32 ±0.05 


0.308 ±0.013 ± 0.036 ±0.006 


B -> /ry 


0.26 ±0.03 


0.24 ±0.04 


0.27 ± 0.011 ±0.032 


B s -> <fo 


0.27 ±0.03 


0.29 ±0.05 




B s -> K * 7 


0.23 ±0.02 


0.20 ±0.04 
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